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In this paper, we study the foUowing question: given a black box performing some unknown 
quantum measurement on a multi-qudit system, how do we test whether this measurement has 
certain property or is far away from having this property. We call this task property testing of 
quantum measurement. We first introduce a metric for quantum measurements, and show that it 
possesses many nice features. Then we show that, with respect to this metric, the following classes 
of measurements can be efficiently tested: 1. the stabilizer measurements, which play a crucial 
role for quantum error correction; 2. the fc-local measurements, i.e. measurements whose outcomes 
depend on a subsystem of at most k qudits; 3. the permutation-invariant measurements, which 
include those measurements used in quantum data compression, state estimation and entanglement 
concentration. In fact, all of them can be tested with query complexity independent of the system's 
dimension. Furthermore, we also present an algorithm that can test any finite set of measurements. 
Finally, we consider the following natural question: given two black-box measurement devices, how 
do we estimate their distance? We give an efficient algorithm for this task, and its query complexity 
is also independent of the system's dimension. As a consequence, we can easily test whether two 
unknown measurements are identical or very different. 
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I. INTRODUCTION 

Quantum measurements are ubiquitous and play a piv- 
otal role in the quantum information science. Many 
tasks require not only the most general quantum mea- 
surements, but also the entangled ones on many parti- 
cles. For example, several quantum data compression 
[ll-Q and communication protocols [6] vi^ork by perform- 
ing some entangled measurement related to the repre- 
sentation theory of symmetric group on a multi-particle 
state. For another instance, certain non-abelian hidden 
subgroup problems are reduced to the problem of effi- 
ciently perfor ming certain joint measurement on multiple 
coset states |7l-[l2|. In order to guarantee the success of 
these tasks in practice, it is crucial to make sure that the 
required measurements are implemented with sufficiently 
high fidelity. 

Now imagine that someone builds a quantum measure- 
ment device and claims that it performs some measure- 
ment that we need. How do we check if this is true? Of 
course, we can use the quantum tomography [l3l - [T9| to 
completely characterize the measurement implemented 
by the device. However, this method is very inefficient, 
especially when the dimension of the system is very large. 
Specifically, a general quantum measurement with k pos- 
sible outcomes on a D-dimensional system is described by 
a collection {Mi, M2, . . . , Mk} of measurement operators 
satisfying the completeness equation 

k 

Y^m}m,^i. (1) 
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Here Mi is a bounded linear operator on the D- 
dimensional Hilbert space and it corresponds to the i-th 
outcome. In order to reconstruct the M^'s, we need to 
determine Q{kD^) parameters. If the system consists of 
n d-dimensional particles (or qudits), then D — and 
hence we must access the device n{k(P"-) times. So this 
approach is practical only for systems of moderate sizes. 

Given the great difficulty of fully characterizing a 
black-box measurement, we need to come up with a bet- 
ter method. Suppose our desired measurement has cer- 
tain property, but the unknown measurement is shown to 
be very different from any measurement possessing this 
property, then wc know that it is surely not what we 
want. Therefore, we can consider devising a test that 
separates the measurements possessing certain property 
from those "far away" from them. Of course, in order 
to claim that two measurements are far away, we need 
to introduce a metric (or distance function) than quan- 
tifies their difference in the first place. We call this task 
property testing of quantum measurements. 

Generally speaking, property testing (20l - [2^ is the task 
of deciding whether an object has certain property or is 
far away from any object possessing this property, given 
the promise that it is one of the two cases. For example, 
given a boolean function / : {0, 1}" — >■ {0, 1} as an ora- 
cle (which receives the input x S {0, 1}" and returns the 
value of f{x)), we may want to determine whether this 
function is linear or far away from any linear functions, by 
querying the oracle as few times as possible. The prop- 
erty testing of classical objects, such as boolean functions 
and graphs, has been extensively studied in computer 
science, and it plays an important role in probabilisti- 
cally checkable proofs (PGP) [23^. Remarkably, many 
properties of boolean functions and graphs are found to 
be testable with very few queries. In fact, sometimes 
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the query complexity is even independent of the object's 
size. The property testing of quantum objects, includ- 

H-Iii], has been 



and accordingly the state of the system after the mea- 
surement becomes 



ing quantum states [2^ and operations 
addressed only recently. It was found that many interest- 
ing classes of quantum states and operations, such as the 
product pure states [l^l and Clifford operations [H, H^l , 
can also be tested with very few copies or queries. 

In this paper, we initiate the study of property testing 
of general multi-qudit measurements. First, we introduce 
a distance function for quantum measurements and show 
that it has many nice features. In particular, if two mea- 
surements are close with respect to this metric, then they 
behave similarly on most input states. Then, we present 
efficient algorithms for testing three interesting classes 
of measurements: 1. the stabilizer measurements, which 
play a crucial role for quantum error correction [30|-|32|; 
2. the fc-local measurements, i.e. measurements whose 
outcomes depend only on a subsystem of at most k qu- 
dits; 3. the permutation-invariant measurements, which 
includes those measurements used in q uantum data com- 
pression [H-Q , state estimation [33l-l35| and entanglement 
concentration (s^ls^j. In fact, all of them can be tested 
with query complexity independent of the system's di- 
mension. Furthermore, we also present an algorithm that 
can test any finite set of measurements. Finally, we con- 
sider the following natural question: given two black-box 
measurement devices, how do we estimate their distance? 
We give an efficient algorithm for this task, and its query 
complexity is also independent of the system's dimension. 
As a corollary, we can easily test whether two unknown 
measurements are identical or very different. 

The remainder of this paper is organized as follows. In 
Sec. II, we introduce a metric for quantum measurements, 
formally describe our problem and also present several 
useful tools for our work. Then, in Sec. Ill, IV and V, 
we study the testing of stabilizer measurements, A:-local 
measurements and permutation-invariant measurements. 
Next, in Sec. VI we give an algorithm that tests any finite 
set of measurements. After that, in Sec, VII we present 
an algorithm for estimating the distance between any two 
measurements. Finally, Sec. VIII concludes this paper. 



II. PRELIMINARIES 
A. A Metric for Quantum Measurements 

Consider a I?-dimensional quantum system. Let Jio 
be its Hilbert space, and let B{'Hd) be the set of bounded 
linear operators on Hd. Then any measurement with 
k possible outcomes on this system can be described 
by M = {Mi,M2,...,Mk} where Mi e BiTio) and 
X^iLi M^Mi = I. If we perform M on the state p, then 
the probability of obtaining outcome i is 



tr(M,pM/) 



(3) 



It is worth noting that each Mi is determined only up to a 
phase. Namely, if we replace Mi by /3iMi for some /3i S C, 
\f5i\ — 1, we are still describing the same measurement. 
One can easily see this from Eqs.© and ([3]). In order 
to conveniently deal with this equivalence relation, we 
make the following definitions. For any A £ B{Hd), let 
[A] = {l3A : /3 e C, = 1}. Furthermore, for any M = 
{M^}l<^<k, let [M] = {{AMJi<,<fe : e C, |A| = !}• 
Then M and N correspond to the same measurement if 
and only if [M] = [N]. 

For some applications, the post-measurement state is 
of little interest. In such cases, the measurement can 
be more conveniently described by a positive- operator 
valued measure{POYM) {Ei, E2, ■ ■ ■ , Ek} where Ei := 
M^Mi. But in this paper we care a lot about the post- 
measurement state. We would say that two measure- 
ments are identical if both the distributions of outcomes 
and the post-measurement states are exactly the same 
for the two measurements performed on arbitrary state. 

We want to be able to compare any two measurements 
that may have different number of outcomes. In order to 
achieve this, we append infinitely many zero operators to 
M — {Afi, M2, . . . , Mk} and rewrite it as {MijigN where 
Mi — for any i > k. In addition, we use \M\ = k to 
denote the number of nonzero M^'s. 

Now let n be the set of all measurements on the D- 
dimensional system. We want to define a metric (or dis- 
tance function) A : x K that should satisfy the 
following natural conditions: for any M — {MijigN, TV — 

1. A(M, A) > 0; 

2. A(M, A) = if and only if [M] = [A]; 

3. A(M, A) = A(A,M); 

4. A(M, L) < A(M, A^) + A(A^, L). 

Furthermore, note that a measurement M — {A/i}igN 
on the system A can be equivalently viewed as a joint 
measurement M®I := {Mi(g)/}igN on the system A®B, 
where B is any ancilla system, and similarly for A^ — 
{AijigN. Since the distance between M and A should 
be independent of whether an ancilla system is appended, 
or what kind of ancilla system is appended, we need 

5. A(M, A) = A(M(g)/, N®I), where / is the identity 
operation on any finite-dimensional Hilbert space. 

Finally, since we want to use A for property testing, it 
needs to be normalized, i.e. 



p{i) = ir{M,pMl), 



(2) 



6. A{M,N) < 1. 
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Our idea is that A(M, N) should reflect the total differ- 
ences between each pair of Afj and Ni. First, we quantify 
the distance between any A,B £ B{'Hd) as: 



A{A,B) 



inf -^\\A'-B' 
A'elA]^B'elB] y/2D 

inf -^\\A- e'^Bllp, 



(4) 



where 



||yl||p:=,/t.r(Atyl) = 



D 

El- 



(5) 



is the Frobenius norm (or Hilbert-Schmidt norm) for 



A = 



(fixing an orthonormal basis for Hd, any op- 



erator A e B{'Hd) is represented by a D x Z? matrix). 
This metric has been used in Ref. for the property test- 
ing of unitary operations. It possesses the following nice 
features: 

(i) A(A,S) >0; 

(ii) A(A, B) = if and only if [A] = [B]- 

(iii) A(A,B)-A(i?,A); 

(iv) A(A,C) < A(AS)+A(B,C); 

(v) A(A, B) = A(A(g) /, B®I), where / is the identity 
operation on any finite-dimensional Hilbert space. 

Furthermore, note that 

A'(^'S) - 7^ ((A A) + {B,B)-2\{A^B)\), (6) 
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where 



{A,B) := tr(AtB) 



(7) 



is the Hilbert-Schmidt inner product. 

Now we define the distance between M — {Af^ligN and 
N = {iVjieN as 



A(M,iV) 



m<Lx{\M\.\N\} 

' J: A2(M„iV,). 

4=1 



(8) 



Properties (i), (ii), (iii) and (v) of A(A, B) immediately 
imply that A{M,N) satisfies conditions 1, 2, 3 and 5. In 
addition, A(M, N) also fulfills condition 6 because, by 
Eq.dH) and the completeness equation. 



A2(M,iV) = EA2(Af„iV0 

-2\{M,,N,)\) 

^ iGN 

< 1. 



(9) 



Moreover, A(M, iV) also satisfies condition 3. To prove 
this, we use property (iv) of A(v4, B) and the following 
lemma: 

Lemma 1. If Qi, bi and ci are non-negative numbers 
such that Ci < Qi + bi for i — 1,2, . . . , k, then 



(10) 



Proof: We need to prove 

E^?<E«? + E^? + 2 

i—l i—1 i—1 

Since < q < + 6^, we have 



(^a^)(^6|). (11) 



E ^ E(°' + ^»)' = E("? + + 



(12) 



So it is sufficient to show 



1=1 



E "^^^^ 



< 



\ .E«?)(E^')' 
\ i=i i=i 



(13) 



which is exactly the Cauchy-Schwarz inequality. 



Plugging a, = A{Mi, N^), h = A(iV,, Li), Ci = A(Mi, Li) 
for 1 < i < niax{|M|, |iV|, |L|} into lemmafl] we obtain 



A(Af, L) < A(A/, N) + A{N, L). 



(14) 



Now we show that A(Af, iV) refiects the average dif- 
ference between the "behaviors" of M and iV on a ran- 
dom input state. Here by "behavior" we mean both 
the distribution of measurement outcomes and the post- 
measurement states. Specifically, let |'0) be a random 
pure state chosen according to the normalized Haar mea- 
sure. If we perform M = {ALijigN or iV = {TVijigN 
on then the unnormalized post-measurement states 
would be {ALi|?/')}igN or {Ni\ijj)}i^^. We quantify the 
total difference between the two sets of states as 

Y.\\Nm-Nm? = Y.mMiM.,+N}N, 

iSN iGN (15) 

-AL/iV,-7V/AL,)|V'). 
Integrating this quantity over and using the fact 

= l (16) 



and the completeness equation, we obtain 
^||Ai,|i/^)-iV,|V)|r = 2 - ^ 



iV,|V)f = 2-^ERe(A//„7V,- 

^ iGN 



(17) 
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By multiplying each Ni by an appropriate phase (with- 
out changing the measurement being described), we can 
make {Mi,Ni) real and non- negative. Then by Eq.®, 

/ V \\Am - N,m^ = 2 - 1- E m^,m\ 

= 2A2(A/,iV). 

(18) 

Thus, if A{M,N) = e is small, then the expectation of 
E.gN - N,\^)\\^ is 2e2, which implies M,\iP) and 

Nilip) are close on average. Furthermore, let piQip)) = 
|lM,|V')f and qi{\i^)) = llA^ilV-)!!^ be the probability of 
obtaining outcome i when performing M and N on {ip) 
respectively. Then by the fact 

\\u~v\\>\\\u\\-\\v\\\ Vu,^;eC^, (19) 

we have 

EiiM#)-iv,iv^)i|2 > Eiimm-imm' 

= 2^-2F(p(|V)),q(|^))), 

(20) 

where 

F{p{m,q{\m--^EVMi^Mm (21) 

is the fidelity of p(|-0)) (Pi(|?A)))ieN and q(|?A)) 
(qi(|-0)))iGN- Taking the expectation of Eq. ipO)) . then we 
know from Eq.dlHl) that F{p, q) is at least 1 - A'^{M, N) 
on average. So if A(M, N) is small, then p and q are also 
close on average. 

Note that by using the Markov inequality 

Pr(|X| >a)< ^1^, Va > 0, (22) 
a 

we can estimate the fraction of "good" input states 
on which M and N behave similarly. Specifically, let 
A(M, N) = e. Then setting X = E.eN ||M,|V')-iV#)|P 
and a = lOe^ in Eq.(l22l) yields 

Pr[^ \\A'm - N,m' > lOe^] < 1 (23) 

So for at least 4/5 fraction of |'0)'s, MIt/i) and A^IV") are 
close, provided that e is small enough. Similarly, setting 
X = l-F(p(|7/')),q(|V'))) anda = lOe^ in Eq.i^ yields 

Pr[l-F(p(|^)),q(|^))) < lOe^] < 1 (24) 

5 

So for at least 4/5 fraction of \tp)'s, p{\i^)) and q(|?A)) are 
close, provided that e is small enough. 

B. Property Testing of Quantum Measurements 

The task of property testing can be formally described 
as follows. Suppose 51 is a class of mathematical objects 



equipped with a metric A. A property is just a subset 
S C ft. For any A G il., ii A G S , then we say that A has 
property S; otherwise, if A{A,S) > e, i.e. A{A, B) > e 
for any B G S, then we say that A is e-far from property 
S. An algorithm e-tests property S with query complex- 
ity (7(|r2|,e) if on any input A £ Q, it makes at most 
(7(|17|,e) queries to A and behaves as follows: 

• HA has property S, then the algorithm accepts A 
with probability at least 2/3; 

• if ^ is e-far from property S, then the algorithm 
accepts A with probability at most 1/3. 

Note that by repeating this algorithm many times and 
choosing the majority answer, we can exponentially de- 
crease the probability of making an erroneous decision. 
So the completeness error 1 — 2/3= 1/3 here and sound- 
ness error 1/3 here can be replaced by arbitrarily small 
constants. 

In this paper, we study the property testing of quan- 
tum measurements on finite-dimensional systems with re- 
spect to the metric A defined by Eq.®. Specifically, 
fix a ZJ-dimensional system, and let be the set of all 
measurements on this system. Suppose 5 C O is the 
property to be tested. Our input is a black box per- 
forming some unknown measurement M G fl. We can 
access it by preparing the D-dimensional system (de- 
noted by A) in some known state pA, applying M on the 
state, and obtaining an outcome i as well as the post- 
measurement state MipaMI (up to a normalization). 
Moreover, we can also introduce an ancilla system B and 
prepare the joint system A® B va some entangled state 
Pab, then apply M on subsystem A of this state, and 
get an outcome i as well as the post-measurement state 
{Mi (g) I)pab{M^ 'S) I) (up to a normalization). An algo- 
rithm e-tests S with query complexity q{D, e) if for any 
input AI G n, it accesses the black box in the above way 
at most q{D, e) times and makes a correct decision about 
whether M e 5 or A{M, S) > e with probability at least 
2/3. Note that we allow the algorithm to extract useful 
information about M from both the outcome statistics 
and post-measurement states. Our goal is to devise such 
a testing algorithm with the minimal query complexity. 
In addition, we also prefer this algorithm to be efficiently 
implementable. Namely, we want its time complexity to 
be polynomial in log(Z3) and 1/e, assuming one query to 
the black box takes a unit time. 



C. Useful Tools 

The following tools will be very useful for our work. 

1. Probability Theory 

Due to the probabilistic nature of quantum measure- 
ments, our testing algorithms heavily depend on analyz- 
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ing the outcome statistics. So the foUowing results from 
probabihty theory will be very helpful. 

Suppose p = (pi)ieN and q — (qi)igN are two proba- 
bility distributions over N. Let 



^(p,q) 



be the variational distance of p and q, and let 



/PiQi 



(25) 



(26) 



be the fidelity of p and q. It is easy to see that 
< A(p, q), F(p, q) < 1. In addition, they satisfy the 
following relation: 



1 - F(p, q) < D{p, q) < ^1-F^p,q). 



(27) 



Meanwhile, the Chernoff-Hoeffding bound [3^, 123] is 
another important and useful result from probability the- 
ory. It gives an exponentially decreasing bounds on tail 
distributions of sums of independent random variables. 
Specifically, let Xi,X2,--.,Xn G [0,1] be independent 
random variables and let S — J27=i-^i- Then for any 
< e < E[S]/n, 



Pr 



> e 



< 2e" 



(28) 



In particular, if the X^'s are independent and identically 
distributed (i.i.d.) as X, then 



Pr 



> e 



< 2e 



(29) 



By choosing n = 0(log(l/i5)/e^) we can make the right- 
hand side smaller than S. Then E[X] can be approxi- 
mated by S/n with precision e and confidence 1 — (5. 

Now imagine that X is a {0, l}-valued random vari- 
able indicating whether a particular outcome occurs in 
an experiment. Namely, X = 1 if a this outcome oc- 
curs, and X = otherwise. Then F,[X] is the probability 
for this outcome in this experiment. By Ea. (|29l) . we can 
approximate this probability with precision e and confi- 
dence 1 — (5 by repeating the experiment 0(log(l/<5)/e^) 
times. 



2. The Choi-Jamiolkowski Isomorphism 

The Choi-Jamiolkowski isomorphism [i^, |4l| states 
that there is a duality between quantum operations and 
quantum states. Specifically, let 



l$1 



(30) 



be the D-dimensional maximally entangled state. For 
any A g B{'Hd), define 



\v{A)) := (A® /)!$+) 



(31) 



where A acts on the first subsystem. Then for any A, B, 



(viAMB)) ^ ^{AB). 



(32) 



Namely, the inner product between |w(A)) and \v{B)) is 
proportional to the Hilbert-Schimdt product of A and B. 
In particular, ii A — B, then we get 



p{A) ■.^\\\v{A))r = ^\\A\\l. 



(33) 



So |t'(^)) is not normalized unless \\A\\f — \/D. We use 
|z;(A)) to denote the normalized version of \v{A)), i.e. 



HA)) 



HA)) 
\v{A))\\ 



(34) 



Now suppose that we perform a measurement M = 
{MijigN on the first subsystem of |$^), then the unnor- 
malized post-measurement states would be {\v{Mi))}i^fi. 
Namely, the outcome i occurs with probability 



p{M,) = \\\v{Au))r = ^\mrF, 

and the corresponding post-measurement state is 

1 



|i)(M,)) = 



HM,)). 



(35) 



(36) 



This property will be crucial for every testing algorithm 
given in this paper. 



3. Pauli Decomposition 



Let 



(37) 



be the Pauli operators. For any x = {xi,X2, ■ ■ ■ ,Xn), 

Z = (^^1, 2^2, ■■ ■ ,Zn) € Z2, let 



Then for any a, b, c, d e Z2 , we have 

[o'a.bO'cd] = [c^aecbed] 



(38) 



(39) 



Here '©' denotes the bitwise addition modulo 2. So there 
exists /3a,b,c,d G C with |/?a,b.c,d| = 1 such that 



O'a.bO'cd — /3a,b,c,dCa©c,b©d- 



(40) 



Moreover, {crx,z}x,zezj form an orthogonal basis for 
Binf") with respect to the Hilbert-Schmidt product. So 
for any A £ B{Hf"') we can write it as 



A^ Mx,z(^)o-x 



(41) 



6 



where 



Note that 



Mx,z(^) := -^{(yyi,^:, A) . 



(42) 



(43) 



Furthermore, by the Choi-Jamioikowski isomorphism, 

(44) 



E A^x,z(^)|«(^x.z)), 



where {\v{ax,z))}x,zeZ'^ form an orthonormal basis for 

(45) 



nf". So 



p{A) = \\\viA))r= E lMx,z(A)p. 



It follows that 

i^(^)> 



^MA) 



KA)) 



^==b(<^x,z)). 



(46) 



So if we measure |t5(A)) in the basis {|i'(crx,z))}x,z6Zj, 
then the outcome (x, z) occurs with probabihty 



q^AA) 



Imx,z(A)|- 

p(A) 



(47) 



Now consider any measurement AI = {M^jigN on an n- 
qubit system. By plugging Mi = Ex.zszj A'x,z(Mi)(Tx,z 
into EiGN M-Mi — I, we obtain 



E E Imx,z(m,)|2 ^ ^ 



(48) 



Furthermore, if we measure \v{Mi)) in the basis 
{|w((Tx,z))}x,zezj, then the outcome (x, z) occurs with 
probabihty 



9x,z(Mj) 



Ia'x,z(m»)|^ 



(49) 



These properties wih be very useful for the testing of the 
stabilizer and fe-local measurements. 

Finally, the above facts can be straightforwardly gen- 
eralized to the qudit case. We mainly need to replace the 
Pauli operators by their higher-dimensional analogues 



J=0 



(50) 



for X, z G Zrf. Here uj = e^'^^/'^ and denotes the addi- 
tion modulo d. The reader can easily work out the rest 
of the details. 



III. TESTING THE STABILIZER 
MEASUREMENTS 

Equipped with the notations introduced above, we be- 
gin with the testing of the stabilizer measurements. This 
kind of measurements ha ve p layed a crucial role in quan- 
tum error correction [30-32]. They are closely related to 
the stabilizer codes ^42], which are defined by choosing 
a set of commuting operators from the Pauli group on n 
qubits. Specifically, let P„ :— {±crx,z : x,z G Z2}. Then 
any operator g G P„ has eigenvalues -1-1 and —1, and 
its (+l)-eigenspace and (-l)-eigenspace are both 2"-"^- 
dimensional. Now if we choose a set of commuting oper- 
ators 31,32, • ■ • G Pn such that -/ ^ (51,32, ■ • • 
then the simultaneous (-l-l)-eigenspace of 3i,S2,---,5fc 
defines the coding space of a stabilizer code. The error- 
correction for this code is quite simple: we simply mea- 
sure each of 31, 32, ■ • ■ , 3fc, and from the outcomes we can 
determine the error, and then we perform the appropriate 
recovery operation. 

A stabilizer measurement is the projective measure- 
ment onto the (±l)-eigenspaces of Ca.b for any a, b G . 
Formally, it is denoted by P(a, b) = {Pi(a, b), P2(a, b)} 
where 



Pi(a,b) 
^'2(a,b) 



O-a.b 



2 ' 



(51) 
(52) 



So M = {MijigN is a stabilizer measurement if and only 
if there exist a, b G such that 



Mo,o(Mi) = Aia,b(A^l) = Mo,o(M2 



-Ma,b(M2 



1 

2' 



Mx,y(M,) = 0, > 3 or (x, y) ^ (0, 0), (a, b). 

. . 

The following lemma says that if the above condition is 
approximately fulfilled, then M = {Mjjjgr} is close to a 
stabilizer measurement. 

Lemma 2. For any < 7,(5 < 1/4 and a, b G Zj, if a 

measurement M = {AfjjigN on n qubits satisfies 



lMo,o(Mi)|,|/ia,b(Mi 



|/iO,o(M2)|,|Aia,b(M2)| G 



4-7,^4+7 



and 



Re(AiS,o(Mi)Aia.b(Afi)) >^-S, 
Re(/iS.o(A'/2)Ma,b(M2))<-i + <5, 



the 



A(M,P(a,b)) < v/87 + 26. 



, (54) 
, (55) 

(56) 
(57) 

(58) 
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Proof: See Appendix 1X1 ■ 

With the help of lemma [U we obtain the following 
result on testing the stabilizer measurements. 

Theorem 1. The stabilizer measurements on n qubits 
can be e-tested with query complexity 0(l/e''). Fur- 
thermore, the testing algorithm can be efficiently imple- 
mented. 

Proof: Given a black box performing some unknown 
measurement M = {Mijigpj on n qubits, we run the 
following test on it: 



Algorithm 1 Testing the stabilizer measurements 



1. Let D = 2", L 



\0.99N] , W = 



20000 



TV 



64 
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2. Prepare L copies of 

(11 ,l2T-','tn)€^2 



\ii,i2, ■ ■ ■ ,in)\ii,i2, ■ ■ ■ ,in)- 



3. Perform M = {MijigN on the first subsystem of each 
copy of If we obtain only outcomes 1 and 2 in 
these measurements, and the fraction of outcome 1 is 

1 , 1 , 

between g ~ g4 ^ 2 64' ^ contmue; otherwise, 

reject M and quit; 

4. Now we should have at least A'^ copies of \v{Mi)) and 
\v{M2)) which are the post-measurement states corre- 
sponding to the outcome 1 and 2 respectively. 

5. Choose T copies of |t;(Mi)) and \v{M2)}, and measure 
each of them in the basis {|i'(o-x,z))}x,zsz5 . If we ob- 
tain only two different outcomes (0, 0) and (a, b) for 
some a, b G Z2 in these measurements, and the fraction 

1 1 

of outcome (0, 0) is between and — \ for 

^ ' ^ 2 64 2 64 

both those measurements on \v{Mi)) and those mea- 
surements on \v{M2)), then continue; otherwise, reject 
M and quit. 

6. Select another W copies of \v{A'Ii)) and \v{A'l2))- For 
each copy, do the following; 

(a) For j = 1, 2, . . . , n, if aj — bj = 0, then set Sj = 
1; otherwise, measure the j-th qubit of the first 
subsystem in the o'aj,bj basis and set Sj to be the 
outcome +1 or —1. 



(b) If the measured state is |i}(Afi)) and Yl 



-1, 



or the measured state is \v{M2)) and Yl ^3 ~ li 

i=l 

then reject M and quit. 
7. Now M has passed all the above tests. Accept M. 



For correctness, we need to prove: 
(1) if M is a stabilizer measurement, then this algorithm 



accepts M with probability at least 2/3; 
(2) on the other hand, if this algorithm accepts M with 
probability at least 1/3, then M is e-close to some sta- 
bilizer measurement. (Taking the contrapositive, we get 
that if M is e-far away from any stabilizer measurement, 
then this algorithm accepts it with probability at most 
1/3.) 

Before proving the two statements, observe that: 

1. In step 3, for i = 1,2, the outcome i should occur 
with probability p{Mi). So, by Eq.d^H) and our 
choice of L, the fraction of outcome i is e^/64-close 
to p{Mi) with probability at least 0.95. 

2. In step 5, for i — 1,2, if the measured state is 
\v{Mi)), then the outcome (0.0) (or (a, b)) should 
occur with probability qo,o{Mi) (or (?a,b(Afi)). So, 
by Eq. (P5)) and our choice of T, the fraction of out- 
come (0,0) (or (a, b)) is e^/64-close to (?o,o(Afi) (or 
qa.,h{Mi)) with probability at least 0.95. 

3. In step 6, for i — 1,2, if the measured state is 
\v{Mi)), then the sj's depend only on the reduced 
state of the first subsystem of \v{Mi)), which is 



Pi 



tr (M/M,) 
Then it is easy to see 



(59) 



pi-[n 5, = i|pi] 



tr(Pi(a,b)pi) 

tr(Pi(a,b)MiMt) 
tr(M/Mi) 



(60) 



Pr[n s,^-l\p2] 



tr(P2(a,b)p2) 

tr(P2(a,b)M2Mj) ^^^^ 
tr(M2^M2) 



If M passes the test with probability at least 1/3, 
then we must have 



Pr[n 



4 ' 
^2 



(62) 
(63) 



because, if otherwise, step 6 would reject M with 
probability at least l-(l-eV4)'*' > l-e"^"'^/'' > 
0.9 by our choice of W. 

Now let us prove statement (1). Suppose M = P(a, b) 
for some a, b e Zj. Since 



p(Pi(a,b))=p(P2(a,b)) = l/2, 



(64) 
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by observation 1 step 3 rejects P(a, b) with probability 
at most 0.1. Next, since 



'7o,o(f'i(a,b)) = ga,b(-Pi(a,b)) = i, 



go,o(P2(a,b)) =qa.b(P2(a,b)) 



2' 



(65) 
(66) 



by observation 2 step 5 rejects /'(a, b) with probabihty 
at most 0.1. Then by observation 3 step 6 never rejects 
P(a, b). So, overaU, algorithm [T] accepts P(a, b) with 
probability at least 0.8. 

Next, we prove statement (2). Suppose M is accepted 
by algorithm [T] with probability at least 1/3. Then we 
claim that: 

(i) M satisfies Eqs.dM]) and ([55]) with 7 = e^/l&; 

(ii) M satisfies Eqs.dSH) and dST]) with 5 6^/4. 
Then, by lemma [21 M is e-close to P(a, b). 

To prove (i), it is sufficient to show 



and 



p{Mi),p{M2) e 

9o,o(Mi),ga,b(A/i) e 
go,o(M2),ga,b(M2) e 

because if they are true, then 



32' 2 



2 1 



e 
32 



1 

32' 2 
32' 2 



32 

£21 

32 



(67) 

(68) 
(69) 



lAix,z(M,)P = p(M,)gx,z(M,) 



C 



1- — ] (I — 

2^ 32) ' V2 ^ 32 

4 16' 4 ^16 
V^ = 1,2, V(x,z) = (0,0),(a,b) 



(70) 



Eq. (j67| holds because, if otherwise, then by observation 

1 step 3 would reject M with probability at least 0.95, 
contradicting our assumption. Similarly, Eas. (|68p and 
([69| also hold because, if otherwise, then by observation 

2 step 5 would reject M with probability at least 0.95, 
also contradicting our assumption. 

Now it only remains to prove (ii). Plugging 



(71) 



X,Z^Z2 

into Eqs. ([5D|) and ([CT]) gives 



Pr[n 



Pr[n 



llPr] = 2 



1 2Re(AiS,o(A^i)Aia,b(Mi)) + Ai 



2p(Mi) 

(72) 

2Re(^S.o(Af2)/ia,b(Af2)) + A2 



2p(M2) 



where 

Aj := X] ^x,z,x®a,zebMx,z(Afi)Mx®a zeb(-^«) 

(x,z)7t(0,0),(a,b) 

(74) 

for i = 1,2. Note that 

|Az| < J2 |Mx,z(Mi)| • |/Xx®a,z®b(Mj)l 

(x,z)#(0,0),(a,b) 

< E IMx,z(A/,OP 

(x,z)#(0,0),(a,b) 

= p{M,) - \fio.o{Mi)\^ - lMa.b(M,)P 

- \2 32I \A 16 
3e2 

(75) 

where in the second last step we use Egs. lpT)) and (|70p . 
Now by Eqs.dH]), dHT]), d?!]) and ([75]), we obtain 



1 e2 

Re(/xS,o(Mi)Aia,b(Afi)) > 4-7- 



Similarly, by Eqs.(IS21), (EB, GSl) and ([75]), we get 



(76) 



Re(AiSo(Af2)Ma,b(M2)) < - 



1 



4' 



(77) 



(73) 



as desired. 

Finally, algorithm [T] has query complexity 0(l/e'*). 
Moreover, besides querying the black box, it requires 
0(n/e'*) quantum operations including: 1. preparing the 
state which is equivalent to preparing the j-th and 

{n + j)-th qubits in the Bell state for j ~ 1, 2, . . . , n; 2. 
measuring a 2n-qubit state in the basis {|w(crx,z))}x,zGZj, 
which is equivalent to measuring the j-th and {n + j)-th 
qubit in the Bell basis for j — 1,2, ... ,n; 3. measuring 
any Pauli operator on a qubit. In addition, the classi- 
cal processing is easy. So algorithm [1] can be efficiently 
implemented. ■ 

The stabilizer code has been generalized to the qu- 
dit case where d is any prime number. So we can also 
consider testing the stabilizer measurements on n qu- 
dits, which are the projective measurements onto the 
d eigenspaces of (Tx,z for any x, z S ZJJ. Algorithm [T] 
can be straightforwardly generalized to test these mea- 
surements. We basically follow the same pattern, except 
that now we have d equally likely outcomes instead of 
two. The reader can easily work out the rest of the de- 
tails. This generalized algorithm still has query complex- 
ity 0(l/e^), and still can be efficiently implemented. 

Finally, note that algorithm 1 does not only test 
whether M is close to a stabilizer measurement, but also 
learns which stabilizer measurement M is closest to. This 
is fortunate, but not generic. Usually learning an object 
is much harder than testing it, as we need to acquire 
more structural information about the object. 
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IV. TESTING THE fc-LOCAL MEASUREMENTS 

Given a measurement on an n-qudit system, it is nat- 
ural ask whether this measurement truly involves all of 
the n qudits. Namely, its outcome might depends only on 
a small subsystem. Specifically, Let [n] := {f , 2, . . . , n}. 
For any T — {ji,j2, ■ • ■ , jm} C [n], we use T to denote 
the subsystem consisting of the ji-th, j2~th, . . . , jm~ 

th 

qubits. We say that a measurement M = {Mi}i^^ is 
k -local if there exists some T C [n] with \T\ = k such 
that 



(78) 



for any i e N. Here Mi acts on the subsystem T and / 
acts on the complementary subsystem T"^. In this section, 
we will study the testing of these fc-local measurements. 

For convenience, we introduce the following notations. 
For any x = {xi,X2, ■ ■ ■ ,Xn),z = {zi, Z2, ■ ■ ■ , Zn) G ZJJ, 
let 

supp(x, z) = {i e [n] : Xi ^ OT Zi ^ 0} (79) 
be the support of (x, z). Then let 

Ft = {(x, z) : X, z e Z2, supp(x, z) C T}. (80) 
Then for any A e B{'Hf"'), we can write it as 

A = /t(A)+5t(A), (81) 

where 

MA) := >^ Mx,z(^)'Tx.z, (82) 



9t{A) 



X! Mx,z(^)0'x.z, 

(x,z)erT 

Mx,z(^)0'x,z- 

(x,z)0rr 



(83) 



While friA) acts non-trivially only on the subsystem T, 
griA) does not. Also, friA) and griA) are orthogonal 
with respect to the Hilbert-Schmidt product. Therefore, 

\\A\\l^\\MA)\\l + \\gTiA)\\l. (84) 

With the above notations, a measurement M = 
{MijigN is fc-local if and only if there exists some T C [n] 
with |r| = fc such that M, = friMi) for any i e N. The 
following lemma says that if this condition is approxi- 
mately fulfilled, then M is close to a fc-local measure- 
ment: 

Lemma 3. For any < S < 1, if a measurement M = 
{Mi}igN on n qudits satisfies 



\\fTiM,)\\l>Dil-S'), 



(85) 



where D = d", then it is S-close to a \T\-local measure- 
ment. 

Proof: See Appendix IB] ■ 



With the help of lemma |3l we have the following result 
on testing the fc-local measurements. 

Theorem 2. The k-local measurements on n qudits can 
he e-tested with query complexity 0(fc log(fc/e)/e^). Fur- 
thermore, the testing algorithm can be efficiently imple- 
mented. 

Proof: Given a black box performing some unknown 
measurement M — {Mi}i^f^ on n qudits, we run the 
following test on it: 

Algorithm 2 Testing the fc-local measurements 



1. Let D = d" , L 



"1200fc 


Kt) 


+ 1 




£2 









2. Prepare L copies of 
1 

7B 



1*5) = 



5Z l«l,«2,. 



, in)- 



3. For j = 1,2, ...,L, perform M — {AfijigN on the 
first subsystem of the j-th copy of then, no 
matter which outcome is obtained, measure the post- 
measurement state in the basis {|u((Tx,z))}x,zez5 ■ Let 
(xj, Zj) be the outcome. 

4. Accept Af if and only if | Uf=i supp(xj, Zj)\ < k. 



For correctness, it is sufficient to show: 

(1) if M is fc-local, then this algorithm accepts M with 
certainty; 

(2) if M is accepted by this algorithm with probability at 
least 1/3, then M is e-close to some fc-local measurement. 

Before prove the two statements, observe that when 
we perform M = {AfijigN on the first subsystem of 
1$^), with probability p{Mi) we obtain outcome i and 
the corresponding post- measurement state is \v{Mi)). If 
this post-measurement state is measured in the basis 
{|w(o'x,z))}x,zeZ"j then the outcome (x, z) occurs with 
probability f/x.zf-^i)- So, in step 3, the probability of 
(xj , Zj ) being some (x, z) is 



(86) 



ex.z(Af) := E P(A^.)9x.z(Mz) 
Hence, the probability of (xj,Zj) being from F^ is 



77T(Af) := E Cx,z(Af) 
(x,z)erT 
= E E lA^x,z(M,)P 

ieN(x,z)err ^°'> 

= ^ E I1/t(a/.)III- 

This holds for any j = 1, 2, . . . , L. 

Now let us prove statement (1). Suppose M is fc-local, 
i.e. A^i's satisfy Eq.(1751) for some T C [n] with |r| = fc. 
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Then for any z G N, /ix,z(-^i) is nonzero only for (x, z) G 
Tx . This implies that in step 3 we have (x^ ,Zj) g Tx 
for any j. So IJ^^j^ supp(xj, z^) C T and M is always 
accepted by algorithm 

Next, we prove statement (2). Suppose M is accepted 
by the algorithm with probability at least 1/3. Then we 
claim that there exists T C [n] with \T\ — k such that 



7?t(M) > 1 - 



Suppose on the contrary that for any such T, rjT{M) < 
1 — e^. Then M must be rejected by the algorithm with 
probability at least 0.9. The reason is as follows. The al- 
gorithm accepts M only if (xi,Zi), (x2,Z2), . . . , (x^jZ^) 
satisfy | IJ^^j^ supp(xj, Zj)| < k. For any such (xj,Zj)'s, 
we can find k numbers 1 < ji < j2 < ■ ■ ■ < jk ^ L such 
that 



U SUpp(Xj,Zj 



U supp(xj 

4=1 



J- 



(89) 



(For each a G IJ^^j^ supp(xj, Zj), we choose a (xj;,Zj. ) 
such that supp(xj. ,ZjJ contains a. Since there are at 
most k elements in Uj^i supp(xj,Zj), k such (xj. ,Zj.)'s 
are sufficient.) Now there can be (^) possibilities for 
ji, j2, ■ ■ ■ , jk, and when the (xj. ,Zj. )'s are determined, 
the total support W := IJ^^-^ supp(xj, Zj) is fixed, and 
then any (xj , Zj ) for j ^ , . . . , must be chosen from 
Tw, with probability is at most 1 — e'^ by assumption. 
Thus, the total probability of getting a valid sequence of 
(xj,Zj)'s is at most (^)(1 - e2)L-/c < j^k^-P(L-k) < g.l 
by our choice of L. 
Now, by Eqs.dHZl) and 



Y.\\}t{M,)\\1>D{1 



(90) 



for |r| = k. Then, by lemma [Sj M is e-close to a fc- local 
measurement. 

Finally, algorithm [5] obviously has query complexity 
0(fclog(fc/e)/e^) as claimed. Moreover, besides querying 
the black box, it requires 0(fc log(fc/e)/e^) quantum op- 
erations including: 1. preparing the state l^"^), which 
is equivalent to preparing the j-th and (n -I- j)-th qu- 
dits in the d-dimensional maximally entangled state for 
J = 1, 2, . . . , n; 2. measuring a 2n-qudit state in the basis 
{|w(CTx,z))}x,zeZ2: which is equivalent to measuring the 
j-th and {n + j)-th qubits in the basis {\v((Tx.z))}x,z^i,ci 
for j = 1, 2, . . . , n. In addition, the classical processing 
is also efficient. Thus, this algorithm can be efficiently 
implemented. ■ 



V. TESTING THE 
PERMUTATION-INVARIANT MEASUREMENTS 

In this section, we consider testing a class of multi- 
qudit measurements possessing certain symmetry with 



respect to the permutations of the qudits. This kind of 
measurements have been widely used in many tasks such 
as quantum data compression [ll-0] , state estimation [ssl - 
[sBj and entanglement concentration [s^, [s^l . Usually we 
perform such measurements on p®" to extract certain in- 
formation from it, or to transform it into a better form, 
where p is some mixed state that we are interested in. For 
example, in quantum data compression, suppose an i.i.d. 
quantum source produces a d-dimensional mixed state p 
each time. We first collect n copies of this state, and 
then compress /O®" into a smaller Hilbert space whose 
dimension is approximately 2""^^''). Here S{p) is the von 
Neumann entropy of p, i.e. if p has the spectral decom- 
position 



then 



x=0 



d-l 



(91) 



(92) 



x=0 



The typical method [J Q is that we perform a projective 
measurement {He, / — H^} on p®", and if the outcome 
corresponds to H^, then we succeed. Here H^ is the pro- 
jection operator onto the so-called "e-typical subspace" . 
This subspace is spanned by all \xi, X2, ■ ■ ■ , a;„)'s satisfy- 
ing 



1 



n \p{xi)p{x2) ■ ■ .p{Xn) 



Sip) 



< €. 



(93) 



Note that this condition depends only on how many XiS 
are 0, 1, . . . , d — \ respectively. Therefore, this e-typical 
subspace is invariant under any permutation of the n 
qudits. In other words, the measurement {He, / — He} 
treats every qudit equally. This property is shared by 
many other measurements, including those used in other 
quantum data compression methods Q-Qi state estima- 
tion (33l - [35| and entanglement concentration [H, [svj . 
Thus, as a preliminary step of testing these measure- 
ments, we can test whether an unknown measurement is 
permutation-invariant in the first place. 

Formally, let S'„ be the symmetric group on n elements. 
For any t G S'„, it can be viewed as a unitary operation 
on n qudits as follows: 



r|0i, 02, ...,(pn) = |0r-i(l), 0r-i(2) 



, . . . , J. j , 



(94) 



where the |0i)'s are arbitrary qudit states. We say an 
operator A G B(T-L®^) is permutation-invariant if 

A^tAt\ VreSn- (95) 
Note that this condition is equivalent to 

A = J2pkEf- (96) 
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for some ^fc G C and Ek € B{Hd)- A measurement M = 
{Mjjigf^ on n qudits is said to be permutation-invariant 
if Mi is permutation-invariant for each i E N. In this 
case, if we first apply a permutation to the qudits, then 
perform M, and at last apply the inverse permutation, 
then the effect is equivalent to directly performing M. 



A. Representation Theory Background 

Our algorithm for testing the permutation-invariant 
measurements depends crucially on the following results 
from the representation theory of symmetric group. Let 
d, n be arbitrary integers. Let 



= {A = (Ai,A2 
> A<j > 0, 



, . . . , Xd) : Ai e Z, Ai > A2 

d 

E Ai = n.} 

i=i 



> 



(97) 

be the set of partitions of n into at most d parts. A 
Young diagram of shape A G Id,n (denoted by ^^(A)) is 
a finite collection of boxes arranged in left-justified rows, 
with the i-th row having A^ boxes for i = 1, 2, . . . , d. Fox 
example, FigH] illustrates the Young diagram of shape 
(5,3,1). 



FIG. 1: The Young diagram of shape (5, 3, 1). 

Each A € Xd.7i (or the Young diagram Y{X)) corre- 
sponds to an irreducible representation Vx of the sym- 
metric group Sn, and also corresponds to an irreducible 
representation W\^d of the general linear group GL{d) of 
d X d invertible matrices. Furthermore, we can calculate 
the dimension of Vx and IF\,d as follows. For a box in 
the i-th row and j-th column of a Young diagram Y(A), 
we define its hook length hook(i,j) to be the number of 
boxes that are in the same row to the right of it or in 
the same column below it plus one. For example. Fig. 
[2] illustrates the hook lengths of each box in the Young 
diagram of shape (5, 3, 1). 



7 


5 


4 


2 


1— 1 


4 


2 


1 






1 











FIG. 2: The hook lengths of each box in the Young diagram 
of shape (5, 3, 1) 



Then the dimension of Vx is given by 



Vx 



n hook(i,j)' 



(98) 



while the dimension of Wx^d is given by 

d + j-i 
WA,<i := 11 7 — w. -N ■ 
{■U)eY(\) hook z, J 



(99) 



For more details about Vx and Wx d, please see Refs.|43l. 

El- ' 

The Schur-Weyl duality |45l - l47l | states that there exists 
an orthonormal basis for 'Hf"', which is called the Schur 
basis and lablled by |A)|aA)|6A)sch f'^'' A G ^d.m < oa < 
wx,d ^ ^, ^ bx < Vx — ^, such that for any r G Sn and 
E e GL{d) we have 

^|A)|aA)|6A)sch = |A)|aA)(l^A(T)|6A))sch, (100) 
i?^"|A)|aA)|6A)sch = \X){WxAE)\ax))\bx)^Am 

In other words, the space "H?" can be decomposed as 



= (WA,d ® VA)sch, 



(102) 



where Wx,d is uiA,d-dimensional and Va is t;A-dimensional, 
such that any t G Sn and E e GL{d) satisfy 

r = (/®yA(r))sch, (103) 
E^n ^ (WxAE) ® I)sch. (104) 

Since Wx.d is a continuous mapping from GL{d) to 
GL{ijux,d), and GL{d) is a dense subset of the d x d ma- 
trices, we can naturally extend Wx.d to all d x d matrices. 
Then for any permutation-invariant A £ B{'Hf"'), if 



let 



Then 



and 



A ■■=J2^'kWxAEk)- 



(105) 



(106) 



^|A)|aA)|6A) 



Sch 



|A)(iA|aA))|6A)s,h> (107) 



^ = (^A<8)/)sch. 

Aeid,„ 



(108) 



On the other hand, it is easy to see that any A £ B{T-Lf"') 
satisfying Eq. (ll08p commutes with any t £ Sn (which 
satisfies Ea. (|103p '). since they act on disjoint subsystems. 
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Hence, A G B{'Hf^) is permutation-invariant if and only 
if it has the block-diagonal form in Eq. (|108p with respect 
to the Schur basis. 

Finally, let Ad^„ := {{\,ax,bx) : A G ld,n,0 < 
a\ < wx.d — 1,0 < b\ < v\ — 1}, and fix a bijec- 
tion h : 'Z''^ ^ ^d.n- Then the standard basis state 
|ii,i2, ...,in) can be also labelled as \X,a\,b\), provided 
that h{ii, 12, ■ ■ ■ ,in) = (A, ax, hx). We use Usch to denote 
the unitary operation that converts the standard basis 
to the Schur basis, and call it the Schur transform [48| . 
Namely, if zj, . . . , i„) {\ax,hx), then 

t/"Sch|«l,«2, • • • = |A,aA,6A)sch- (109) 

Its inverse is called the inverse Schur transform. 

Both C/sch and C/g^j^ can be (approximately) implemented 
with time complexity polynomial in n, d and log(l/e), 
where e is the accuracy parameter [isj . 

B. Our Result on Testing the 
Permutation-Invariant Measurements 



Lemma 4. For anj/ Q < 5 < \, if a measurement M = 
{A/j}jgN on n qudits satisfies 

j2\\m%>Dii-s'), (117) 

where D = d" , then M is 5 -close to a permutation- 
invariant measurement. 

Proof: See Appendix [Cl ■ 

With the help of lemma SI we obtain the following re- 
sult on testing the permutation-invariant measurements. 

Theorem 3. The permutation-invariant measurements 
on n qudits can be e-tested with query complexity 0(1/ e^). 
Furthermore, the testing algorithm can be efficiently im- 
plemented. 

Proof: Given a black box performing some unknown 
measurement AI ~ \Ali\i&i on n qudits, we run the 
following test on it: 



It will be convenient to establish the following nota- 



tions. For any A e 8(1-1^"), it can be written as 
E (|A)(A'|(8AA.A')sch 



(110) 



for some Ax^y's. We are particularly interested in the 
diagonal term Ax,x G B{'Wx,d ^ Va). Let {g\,j}o<j<vl-i 
be the generalized Pauli operators on Va such that gx^ — 
I. Then we can expand Ax.x as 



Ax,x = iA «) / + E ^A,i «) gxj 

3^0 



(111) 



for some Ax and ^aj's. So A can be written as the sum 
of three terms: 



where 



A^{A + A + A)sch, 
A := Aa®/, 

A := |A)(A'|®Aa,a'. 



(112) 

(113) 
(114) 
(115) 



While A is permutation- invariant, A and A are not. Fur- 
thermore, A, A and A are mutually orthogonal with re- 
spect to the Hilbert-Schmidt product. Thus, 



\A\\l = \\A\\l + \\A\\% + \\A\\l. 



(116) 



With these notations, a measurement M — {Mi}i^m 
is permutation-invariant if and only if Mi — (Mi)gch for 
any i G N. The next lemma shows that if this is approx- 
imately true, then M is close to a permutation-invariant 
measurement. 



Algorithm 3 Testing the permutation-invariant mea- 
surements 

1. Let D = d" , L = 

2. Repeat the following test L times: 

(a) Prepare 



(n,...,»„)ez5 



il, 12, . . . , in)\ii,i2, 



,in). 



(b) Perform the Schur transform on the first subsys- 
tem of 

(c) Perform M on the first subsystem. 

(d) No matter which measurement outcome occurs in 
step (c), perform the inverse Schur transform on 
the first subsystem. 

(e) Re-label the standard basis of both subsystems as 
|A,aA,feA)'s via the bijection /i : ZJJ — )■ Ad,„. 

(f) Measure the two |A) registers. If the two outcomes 
are not equal, then reject M and quit. 

(g) Measure the two \b\} registers in the basis 
{|^(ffA,j))}o<j<^2 -1- If the outcome corresponds 
to j 7^ 0, then reject M and quit. 

3. Now AI has passed all the above tests. Accept M. 



For correctness, we claim that this algorithm accepts 
M with probability (-^ EieN ll-^^j|lF)^- Assuming this is 
true, then: 

(1) if M is permutation-invariant, then Mi — (Afi)sch for 
any i G N. Then by E^gn mJm^ = /, we obtain 



(118) 
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which imphes that algorithm[3] accepts M with certainty; 
(2) on the other hand, if M is accepted by algorithm [3] 
with probabihty at least 1/3, then by our choice of L we 
must have 

Y,\m\l>D{l-e^) (119) 

(because if otherwise, then the probability of M being 
accepted is at most (1 — e^)^ < e^"^^^ < 0.1, contradicting 
our assumption). Then, by lemma IH M is e-close to a 
permutation- invariant measurement . 

Now we prove the above claim by showing that 
M passes each iteration of step 2 with probability 

^-gpj ||Mi|||,. Assuming the outcome at step 2.(c) is 
«, the unnormalized state after step 2.(e) is given by 

[(C/L^/,C/sch)®/]|$i) 
= [{M, + M, + Mi) ® 

= ^ E {{]^'h + M, + M,)\\axM))\\ax,hx), 

(120) 

where 

m\X)\ax)\bx) = \X}{{m)x\ax})\bx}, (121) 
m\X}\ax}\bx) = E \X)m)x,j\ax)){9x.j\bx}), (122) 

m\X)\ax)\bx) = EJX'){mx',x\ax)\bx)). (123) 

Note that Mi changes |A) and M changes |6a). Only 
Mi leaves both |A) and |5a) intact. In step 2.(f), we do 
not reject M only if the measurements on the two |A) 
registers yield the same outcome. In this case, assuming 
the outcome is A, then the unnormalized state of the rest 
four registers becomes 

4= E [m)x\ax))\bx)\ax)\bx) 

VJJ ax,bx 

+ E{mxM)i9xM)\ax)\bx)] 




(124) 



+ Eimx>x))M\H9x,M 

where in the first step we have switched the order of the 
first |6a) and the second \ax)- In this case, the mea- 
surement in step 2.(g) produces the outcome j = with 
probability 

j:\\i:m)x\ax))\ax)r = ^ii(m.)a|ii. (125) 



Summing this probability over all i's and A's, we know 
that M passes one iteration of step 2 with probability 

EE||K^^)aIIf = ;^EII^^^II^' (126) 
as claimed. 

Finally, this algorithm obviously has query complex- 
ity 0(l/e^). Furthermore, besides querying the black 
box, it needs poly(n, 1/e) number of quantum operations 
including: (l)preparing the state |$^), which is equiva- 
lent to preparing the j-th and (n + j)-th qudits in the 
d-dimensional maximally entangled state; (2)quantum 
Schur transform and its inverse, which can be efficiently 
implemented; (3)measuring the two \bx) registers in the 
basis {\v{gx,j))}o<j<v\-i, which can be accomplished ef- 
ficiently with the circuit in FigH) 



FIG. 3: The circuit for measuring the two |6a) reg- 
isters in the basis {|«((/A,i))}o<j<^2 _j where the (/aj's 
are the wa- dimensional Pauli operators. Here CXv^ is 
the HA-dimensional generalization of CNOT: CXv^\a)\b) = 
\a)\b + a(mod v\)), for a, 6 = 0, 1, . . . , ua — 1, and Fv^ is the 
UA-dimensional quantum Fourier transform. Since t;A < d", 
both CXv^ and F^^ can be efficiently implemented. 

In addition, the classical processing is also easy. So 
algorithm [3] can be efficiently implemented. ■ 

VI. TESTING ANY FINITE SET OF 
MEASUREMENTS 

So far we have studied testing three particular classes 
of quantum measurements. In this section, we present 
a general result about property testing of quantum 
measurements. Specifically, we give an algorithm that 
can test any finite set of measurements on any finite- 
dimensional system. 

Before stating our result, it is helpful to consider 
the following question. Suppose M — {Mi}i^fi and 
N = {A^ijigN are two measurements on a _D-dimensional 
system. If a black box performs one of them, how do we 
know which one it really implements? We solve this prob- 
lem by repeating the following experiment sufficiently 
many times: each time we prepare a copy of l*!"^) and 
apply the unknown measurement to its first subsystem, 
obtaining an outcome and post-measurement state. If 
the unknown measurement is M (or A^), then the out- 
come i occurs with probability p{Mi) (or p{Ni)) and the 
corresponding post-measurement state would be \v{Mi)) 
(or \v{Ni))). We claim that as long as A(M, N) is large, 
we do not need to repeat this experiment many times 
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to successfully identify the unknown measurement. The 
reason is as follows. 

Let p(M) := (p(M,))»eN and p(iV) ip{N,)),eN be 
the distributions of outcomes for M and N respectively. 
Consider their variational distance 

A(p(Af),p(7V)) ^ ^J2\p{M,) - p{N,)\. (127) 

By Eqs.®, (HZ]), and the Cauchy-Schwarz inequal- 
ity, 

A^M,N) = 1- E \{v{M,MN,))\ 

> i-EI!KM.)>ll-IIKA^.))ll 

= 1 - E VpiMMNT) (128) 
= l-F{p{M),p{N)) 

> ^D^{p{M),p{N)). 

Thus 

A(M, N) > ^D{p{M), p{N)). (129) 

This means that if D{p(M),p{N)) is large, then 
A{M,N) is also large. In this case, we can simply dis- 
tinguish M and N from the outcome statistics. 

But what if A{M,N) is large but L»(p(M), p(A^)) is 
small? In this case, we need to use the post-measurement 
states. Specifically, suppose \M\, \N\ < k. We repeat the 
above experiment L times, and suppose the outcome i 
occurs Li times for i = 1,2, ... ,k. Note that with high 
probability we have Li « p{Mi)L « p{Ni)L for each i. 
Then let L = (Li, L2, . . . , ife) and 

k 

Xl(M) :^ (g)(|5(M,))®^0, (130) 

k 

XL(iV):^(g)(|5(iV,))®'^')- (131) 

The following lemma says that if A{M, N) is large, then 
with high probability xl (M) and xl (N) have small over- 
lap for some small L, and thus can be easily distin- 
guished. 

Lemma 5. Let L = {Li, L2, ■ . ■ , Lk) be a set of non- 
negative integers with Ei=i ~ ^- Suppose M ~ 
{Ml, M2, . . . , Mk} and N ^ {Ni,N2, . . . , iVfc} are two 
measurements such that A{M,N) > S for some < S < 
1. If for any i such that Li > O.lS^L/k, we have 

p{M,),p{N,)> (1-0.15^)^, (132) 

then 

|(XL(Af)|xL(iV))|<(l-0.6<52)^ (133) 



Proof: See Appendix ID] ■ 



Now let us return to the problem of testing any fi- 
nite set of measurements on a finite-dimensional system. 
Suppose S = {M(i), M(2), . . . , m(")} is the set of mea- 
surements to be tested, where M*^*^ = {Mj'''}jgp}- Given 
an unknown measurement M, we want to know whether 
AI is in S or far away from S. Our basic idea is to repeat 
the above experiment sufficiently many times. If M is 
from S, then with high probability the outcome statis- 
tics is well-behaved and the post-measurement states lie 
in certain subspace. On the other hand, if M is far from 
any M^^\ then we divide S into two subsets: for one sub- 
set, the outcome statistics are quite different for M and 
any M*-*^ in it, so we can easily distinguish M from this 
subset; for the other subset, the outcome statistics are 
similar for M and any M*^*) in it, so we must utilize the 
post-measurement states. Lemma[5]says that XL(Af) and 
XL (Af *■*■*) are almost orthogonal for some small L. Then 
by the following lemma, xl(M) has a small projection 
onto the subspace spanned by these xl(-W^*^)'s. So by a 
projective measurement, we can get to know whether M 
is in S or far away from S. 



Lemma 6. Suppose l^i), |02), . . . , |(/>m), IV") G Hd 
(where D is arbitrary) satisfy \{ip\(l)i)\ < l/(5m) for any 
i and \ {4>i\4>j)\ < l/(5m) for any i ^ j. Let 11 be the pro- 
jection operator onto the subspace spanned by \4>i) , \4'2) , 
|0„,). Then (^|n|^) <0.1. 



Proof: See Appendix IE] ■ 



Now we formally state our result on testing any finite 
set of measurements on any finite-dimensional system. 



Theorem 4. Suppose S = {M^^\M^'^\. . . , M^™)} is a 
set of measurements on a D-dimensional system, where 
AfW = {M^\m^\...,mI^^} has at most k possible 
outcomes for any i, and A(MW,M(J')) > 7 for any 
* 7^ J; for some < 7 < 1. Then S can be e-tested 
with query complexity 0(max{fc^ log(fc)/a^, log(r7i)/a^}) 
where a = min{e,7}. 



Proof: Given a black box performing some unknown 
measurement M — {AfjjjgN, we run the following test 
on it: 
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Algorithm 4 Testing S = {M^^), M^^), . . . , M^™)} 



1. Let L — max 



5000fc^ ln(20fc) 



21n(5m) 



D-l 

„ E 1^)1^). 

D 1=0 



2. Prepare L copies of j^^) 

3. Perform M = {MijigN on the first subsystem of each 
copy of 1$^). Let Lj be the number of occurrences of 
outcome j for any j £ N, and let L — (Li, L2, . . . , Lfc). 

4. If Lj > for some j > k, then reject M and quit; 
otherwise, define T C 5 as follows: M^*' € T if and 
only if for any j such that Lj > O.la^L/k, we have 
p(Mj'') > (1 - 0.la'^)Lj/L. If r = 0, then reject M 
and quit. 

5. Now suppose r = {M('i\m('2),...,m(**'} for some 
1 < t < m. Let n be the projection operator 
onto the subspace spanned by Xl(-^^*^^)i 

. . . , Perform the projective measurement 

{II,/ — n} on the state K the outcome corre- 

sponds to n, then accept M; otherwise, reject M. 



Before prove the correctness of this algorithm, observe 
that, in step 3, we should obtain outcome j with proba- 
bility p{Mj), for any j. Then by Eq. (|29| and our choice 
of L, 



Pr 



■p{M,) 



< 



lOOfc 



> 1 



1 



(134) 



for any j. Thus, with probability at least 0.9, for any j 
satisfying Lj/L > O.la^/fc, we have 



lOOfc - ^ ' L 



(135) 



Now suppose that M — M^*) for some i. By the above 
observation, M*^*^ is in the set T with probability at least 
0.9. Then, xl{M) = xl(M(*)) is obviously in the sub- 
space spanned by the XL(Af(*J^)'s and hence the measure- 
ment in step 5 always produces the outcome correspond- 
ing to n. Therefore, M = M^^^ is accepted by algorithm 
IHwith probability at least 0.9. 

On the other hand, suppose AI is accepted by algo- 
rithm 2] with probability at least 1/3. We need to prove 
that M is e-close to some M*^*^. Suppose on the contrary 
that this is not true. Then, still by the above observa- 
tion, with probability at least 0.9, for any j satisfying 
Lj/L > O.laVfc, we have p{Mj) > (1 - 0.1a'^)Lj/L. Fur- 
thermore, by the definition of T, any M^^'^ € T also 
satisfies p{M^/''>) > (1 - Q.la'^)Lj/L for such j's. Then 
by lemma O and our choice of L, we obtain 

|(xl(M)|xl(M('')))| < (l-0.6a2)^ 

^ g-0.6a^L 
1 



(136) 



< 



5m ' 



for any I, and similarly 

|(xl(M(^'))|xl(M(v)))| < 



1 

5m ' 



for any I 7^ I' . Thus, by lemma[6l in step 5 we obtain the 
outcome corresponding to 11 with probability at most 
0.1. Overall, M can be accepted by algorithm |4] with 
probability at most 0.2, contradicting our assumption. 

Finally, the query complexity of algorithm |4] is 
0(max{fc^ log(fc)/a^, log(m)/a^}) as claimed. In partic- 
ular, if the system consists of n qudits, then the query 
complexity is polynomial (in n) as long as k and 1/a are 
at most polynomial and m is at most exponential (in n). 
■ 

In general, algorithm 2] may be not efficiently imple- 
mentable. But from an information-theoretic point of 
view, theorem 2] shows that we do not need to query the 
black box too many times to know whether it belongs to 
some finite set or is far away from this set, as long as the 
set is not too large (i.e. m is not too large) , the mea- 
surements in the set are well-separated (i.e. 7 is not too 
small), and any measurement in the set does not have 
too many possible outcomes (i.e. k is not too large). 

Finally, let us apply theorem |4] to the set of sta- 
bilizer measurements. This set contains 4" elements. 
The distance between any two of them is 1/2, i.e. 
A(P(a,b),P(c,d)) = 1/2 for any (a, b) ^ (c,d). Fur- 
thermore, each stabilizer measurement has only 2 pos- 
sible outcomes. Thus, by theorem 01 the stabilizer 
measurements can be e-tested with query complexity 
0(max{l/e^, n/e^}). This query complexity is only 
slightly worse than the 0{l/e'^) given by theorem[TJ This 
example shows that although algorithm U uses only cer- 
tain distance information about the measurements to be 
tested, it still can achieve a quite good efficiency. 



VII. ESTIMATING THE DISTANCE OF TWO 
UNKNOWN MEASUREMENTS 

In the previous sections, we have studied testing sev- 
eral properties of a single measurement. In this section, 
we consider a different scenario. Suppose we are given 
two black-box measurement devices, how do we know 
whether they perform the same measurement? And if 
not, how do we estimate their distance? We will give an 
efficient algorithm for this task. Surprisingly, its query 
complexity does not depend on the dimension of the sys- 
tem, but depends only on the proximity parameter and 
the number of possible outcomes for the two measure- 
ments. 

Our basic idea is as follows. Suppose M — {Mi}i<ci<k 
and N = {Ni}i<ci<k are two measurements with at most 
k possible outcomes. Since 

k 

A\M,N) = l-Y. MM,)\v(N,))\. (138) 

i=l 

So we can estimate (M, N) by estimating each 
\{v{Mi)\v{N,))\. Note that 



(137) 



\{v{M,)\v{Ni))\ = ^p{M,)p{Ni)\{v{Mi)\v{N,))\. 



(139) 
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Thus, we can estimate \{v{Mi)\v{Ni))\ by estimating 
p{Mi), p{Ni) and \{i{M,)\v{Ni))\. Recall that p{Mi) 
(or p{Ni)) is the probabihty of obtaining outcome i 
when we perform M (or N) on the first subsystem 
of and \v{Mi)) (or \v{Ni))) is the corresponding 

post-measurement state. So we repeat this experiment 
many times. From the outcome statistics, we can es- 
timate p{Mi) and p{Ni) with good precision. As to 
\(v(Mi)\v{Nj))\, we can estimate it by performing the 
swap test [431 o^i the states \v{Mi)) and \v{Ni)). In fact, 
the swap test can be used to estimate the overlap between 
any two pure states, as stated by the following lemma. 

Lemma 7. Given 0(log(l/i5)/e'') copies of \<f)) and \^), 
we can estimate \{4>\'ip)\ with precision e and confidence 
1 — 5. Furthermore, the estimation algorithm can he ef- 
ficiently implemented. 

Proof: Given any |0) and \ip), the swap test is the 
standard procedure to estimate It works as 

follows: first, we prepare a qubit in the state |-|-) — 
1|^/2{\Q) + |1)); then, we perform a controlled-swap gate 
on 10) and It/;), using |-|-) as the control qubit (a swap 
gate is the operation |0)|V') — ^ IV') 10)); finally, we apply 
a Hadamard gate on the control qubit and measure it 
in the standard basis. The circuit for this procedure is 
illustrated in Fig. IH 





H 






H 




















Swap 





FIG. 4: The swap test 

It is easy to see that the final state before the measure- 
ment is 

\\mm) + + \\mm) - \m))- im 

Thus, the measurement on the first qubit yields outcome 
with the probability (1 -I- 1(01-0) p)/2. 

Now consider the following estimation algorithm: 



Algorithm 5 Estimating K^IV")! 

1. Assume we have L = [21n(2/(5)/e''] copies of 10) and 
|'0}. We run the swap test on each pair of |0) and 
Let po be the fraction of outcome among the out- 
comes obtained. Then return ^2po — 1 as the estimate 
of I (01^)1 . 



By Eq. (l29)) and our choice of L, po will be e^/2-closc 
to (1 + |(0|i/))p)/2 with probability at least 1-5. This 
implies that 2pQ — l is e^- close to |(0|'0)p with probability 
at least 1 — 5, and hence \/2po — 1 is e-close to |(0|'0) | with 



probability at least 1 — 5 (note that \a — &p < — 
for any a, 6 G [0, 1].) 

The swap test requires three kinds of quantum opera- 
tions including the Hadamard gate, controlled-swap gate 
and the measurement in the standard basis. If |0) and 
10;) are n-qudit states, then the controlled-swap gate can 
be implemented with 0{n) basic quantum gates. Algo- 
rithm [S] repeats the swap test 0(log(l/(5)/e'*) times. So 
it can be efficiently implemented. ■ 

So, if we have sufficiently many copies of \v{Mi)) and 
|u(iVi)), then we can estimate their overlap with good 
precision. But what if p{Mi) or p{Ni) is so small that the 
outcome i almost never occurs, unless we repeat the ex- 
periment extremely many times? Note that for such i's, 
|jMi||F or llA^illF must be small. The next lemma shows 
that the total contribution of such \{v{Mi)\v{Ni))\'s to 
A^(M, N) is small, and hence we can safely ignore them. 

Lemma 8. Suppose M = {Mi, M2, . . . , Mfc} and N = 
{Ni, N2, ■ . . , Nk} are two measurements. Let Xs — {i : 
p{Mi) <5 or p{Ni) < 5} for any < 5 < 1 . Then 

Y,MM^)\v{N,))\<2V5k. (141) 

Proof: Let A ^ {i : p{M,) < 5} and B = {i : p{N^) < 
5}. Then Is — AU B. So it is sufficient to show 

J2MM,)\vm)\<^, (142) 

jGA 

J2\{viA'hMN,))\<V6k,. (143) 

ieB 

To prove Eq. (|142p . we use the Cauchy-Schwarz inequality 
and the definition of A, 

<V5Y. VpW) 

ieA 

< V5k, 

(144) 

where the last step follows from ^ p{Ni) < 1 and \A\ < 

ieA 

k. Ea. (|143p can be proved similarly. 

■ 

Now we are ready to formally describe our result. 

Theorem 5. Suppose M = {Mi, M2, . . . , M^} and N = 

{A^i, A^2, ■ • ■ T Nk} are two measurements with at most k 
possible outcomes on a D-dimensional system. Then we 
can estimate their distance (i.e. A{M,N)) with preci- 
sion e and confidence 0.8 by using them 0{k^ log(fc)/e^^) 
times. Furthermore, the estimation algorithm can be ef- 
ficiently implemented provided that k is polylogarithmic 
in D. 
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Proof: We run the following algorithm on M and N: 



Using these facts and lemma [U we obtain 



Algorithm 6 Estimating A{M,N) 



1. Let L = 



50000fc^ ln(40fe) 



T = 



(- 

\ 16k 



36fc2 



1 

2. Prepare L copies of |<1?^) = — r= Perform 

VD i=o 

M on the first subsystem of each copy. Let be the 
fraction of outcome i among the L outcomes obtained, 

forz = l,2,...,fc. LetA = {^:a.>^-3^}. 

3. Prepare another L copies of |$^). Perform A'^ on the 
first subsystem of each copy. Let hi be the fraction 
of outcome i among the L outcomes obtained, for i = 



1,2, ...,/fc. Let B = {i : &i > 



r}. 



16fc 36fc2 

4. For any i £ Af] B , we have at least T copies of \v{Mi)) 
and \v{Ni)). Then we run algorithm [5] on them, and 
get an estimate of \ {v{Mi)\v(Ni))\. 



5. Return A 



A(M,iV). 



X] V<Hbi\i as an estimate of 

:eAnB 



Now we prove the correctness of this algorithm. Note 
that in step 2, by Eg.ip^ and out choice of L, Oj 
is e^/(36fc^)-close to p{Mi) with probability at least 
1 — l/(20fc) for each i. Therefore, with probability at 
least 0.95, every i satisfying p{Mi) > e'^/{16k) is in A. 
Similarly, in step 3, with probability at least 0.95, ev- 
ery i satisfying p(iVj) > e^/{16k) is in B. Furthermore, 
by the proof of lemma [7] and our choice of T, for each 
i G AnB, the produced by algorithm [5] is e^/(6fc)-close 
to \{v{Mi)\v{Ni))\ with probability at least 1 - l/(10fc). 
Thus, by a union bound, with probability at least 0.8, we 
simultaneously have: 



|A2 - A2(M,iV)| 
- I( E \{viM,MN,))\ - Z v^A,) 

ieAnB ieAPiB 

+ E \{v{M,)HN,))\\ 

i^AnB 

< E \\{viM,)\v{N,))\-V^A,\ + '^ 
ieAnB ^ 

< — + — 



where in the third step we have used 
\\{v{M,)\viN,))\ - V^iM 



(145) 



+ \MM,MN,))\-X4 

e2 £2 g2 
< \ \ 

6^ 6k 6k 

= Vie AnB. 

2k 



(146) 



(Note that for any < ai, a2, /32, 7i, 72 < 1, we have 
|ai/3i7i - a2/5272| < I"! - "2! + |A - /32| + I71 - 72!)- 
Thus, 



|A- A(Af,A^)| < 



(147) 



as desired. 

Finally, this algorithm obviously has query complex- 
ity 0(A:^ log(fc)/£^2)^ Moreover, besides querying the 
black boxes, it requires two kinds of quantum operations: 
1. preparing 0(A:^ log(fc)/£-^2) copies of |$^); 2. running 
algorithm[5]0(fc) times. |<I>^) can be easily prepared, and 
algorithm[5]is also efficient. In addition, the classical pro- 
cessing is also efficient. So as long as k is polylogarithmic 
in D, algorithm [5] can be efficiently implemented. ■ 



1. \a,~p{M,)\ < 



36k^ 



As a corollary of theorem [SJ we can easily determine 
whether two unknown measurements are the same or 
quite different. 



2. \h-p{N.)\<^,W, 



3. p(M,)< — ,Vi^A; 



Corollary 1. Suppose M = {A/i, M2, . . . , Af^} and 
N = {Ni, N2, . . . , Nk} are two measurements with k pos- 
sible outcomes on a D-dimensional system. Assuming 
that they are either the same or e-far away, we can know 
which case it is with probability at least 0.8 by using them 
0(fc^ log(A:)/£^2) times. Furthermore, the testing algo- 
rithm can be efficiently implemented provided that k is 
polylogarithmic in D. 



4. p{N,)<—,yi^B; 



5. |A,-|({}(A/,)|z)(iV,))ll < g^,Vze Ans. 



Proof: We simply use algorithm |6] to estimate A(Af, N) 
with precision £/2. If it is smaller than £/2, then we 
conclude that M and N are the same; otherwise, we con- 
clude that M and N are £-far away. By theorem [SJ we 
succeed with probability at least 0.8. In addition, the 
query complexity and time complexity are as claimed. ■ 
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VIII. CONCLUSION 

To summarize, we have introduced a metric A for 
quantum measurements on finite-dimensional systems. 
This metric indicates the average difference between 
the behaviors of two measurements on a random input 
state. Then we show that, with respect to this met- 
ric, the stabihzer measurements, fc-local measurements 
and permutation-invariant measurements can be aU effi- 
ciently tested with query complexity independent of the 
system's dimension. Moreover, we also present an algo- 
rithm for testing any finite set of measurements. Finally, 
we give an efficient algorithm for estimating the distance 
of two unknown measurements, and its query complex- 
ity is also independent of the system's dimension. As 
a consequence, we can easily test whether two unknown 
measurements are identical or quite different. 

It is worth noting that entanglement plays a crucial 
role in all of our testing algorithms. Namely, we need to 
prepare the maximally entangled state |$^), then "im- 
print" the measurement operator Mi on it, obtaining the 
state \v{M{)) = {Mi ® /)|<i>^), and finally extract in- 
formation about Mi from this state. It seems that by 
utilizing entanglement, we can somehow gain a better 
understanding of the global property of a measurement. 
We suspect that the three classes of measurements con- 
sidered in this paper cannot be tested with the same ef- 
ficiency if entanglement is not allowed. 

Finally, we would like to point out several directions 
that deserve further investigation: 

First, in this paper we have focused on giving upper 
hounds for the query complexity of testing quantum mea- 
surements. It remains open to prove lower bounds for 
the query complexity of the same task. In particular, it 
would be interesting to know whether our testing algo- 
rithms are optimal and, if they are not, how to improve 
them. 

Second, in this paper we have assumed that we have 
access to both the measurement outcome and post- 
measurement state. For some applications, however, 
the post-measurement state is inaccessible and people 
are merely interested the outcome. In such cases, the 
measurement can be more conveniently described by a 
POVM {Ei,E2,...,Ek} where Si > and Y!1=i E^ = I- 
It is worth studying the testing of POVMs. Note that in 
this new problem, we can extract information only from 
the outcome statistics. So it seems to be harder than the 
one considered in this paper, and requires quite different 
techniques. 

Third, it would be also interesting to study the test- 
ing of quantum channels. Quantum measurements and 
quantum channels are closely related. Namely, if we per- 
form a measurement {Mi, M2, . . . , Mk} on a state p but 
the outcome is unknown, then this process can be equiv- 
alently characterized as a quantum channel transforming 
P into MipMj. On the contrary to testing POVMs, 

in testing quantum channels we can obtain information 
only from the mixed output states. So it also seems to 



be more difficult than the one considered in this paper, 
and needs new techniques. 

At last, the query complexity of property testing de- 
pends crucially on the metric used. Here we have consid- 
ered A which is a very natural and reasonable distance 
function for quantum measurements. It would be in- 
teresting to study the testing of quantum measurements 
with respect to other reasonable metrics. 
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Appendix A: Proof of Lemma [5] 

Without loss of generality, we can assume that 
Mo,o(Mi) and /io,o(Af2) are both real and non- negative 
(if not, we can multiply Mi or M2 by appropriate phases 
to make this true). Then by Egs. fM)) and ([55)1 we can as- 
sume iJ,o,o{Mi) = ai, /^o,o(M2) = a2, Pa,h{Mi) = /^le'^i 
and Aia,b(M2) = /32e'^= where 



ai,a2,/3i,^2 e 



4-7,^4+7 



(Al) 



and 6*1, 6*2 G [0,27r). Plugging this into Eqs.dSll) and ^ 
we obtain 

cos(^?i)> 1^1^ > 1-45-47, (A2) 

cos(g2) < - ~ ^ < -1 + 4(5 + 47. (A3) 
1/4 + 7 

Meanwhile, combining Eqs.(|3Hl), (EH) and yields 

E E lMx,z(M,)|2-t- ^ IMx,z(M,)|2 

»=1,2 (x,z)7t(0,0),(a,b) i>3x,zGZ5 
< 47. 

(A4) 

Now consider the distance between Mi and Pi(a, b) for 
different i's: 

A2(Mi,Pi(a,b)) < -L||Mi-Pi(a,b)||| 



ai 



Pie'- 



E 



lMx.z(Mi)|2, 



(x,z)#(0,0),(a,b) 



A2(M2,F2(a,b)) < --||M2-P2(a,b)||| 



(A5) 



a2 



1 



E 



IMx.z(M2)|2, 



(x,z)#(0,0),(a,b) 



(A6) 
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A2(M„P,(a,b)) = ^\m\% 

^ o E lMx,z(AfOP, Vz>3. 
^ x.zezj 

(A7) 

Here L» = 2". Taking the sum of Eqs.dM]), (TM)) and 
[7| for alH > 3 and using Ea. (|A4l) . we obtain 



A^(M,P(a,b)) < - 



1 



1 




1 


2 


1 




1 


2 


Q;2 — 


2 


H 


^2 


^2 



+27 

Now by Ea. (UT|) we have 





1 


2 


2 1 


ai - 


2 


< 


"1-4 




1 


2 


2 1 


0L2 - 


2 


< 


a, 

2 4 



<7, 
<7- 



(A8) 

(A9) 
(AlO) 



Moreover, by Eqs. (IXTt - (|Mt . we have 
1 2 



1 



/3f + T-/?i cos(0i) 



/32e" 



(1 - 47 - 4(5) 
< 57 + 2(5, 



= /3| + 7 + /?2 cos(02) 



(All) 



(1 - 47 - 4(5) 
< 57 + 2(5. 



(A12) 



(A13) 



Plugging Eqs.(R9l)- (|ll2l) into Eq.jM]) yields 
A2(Af,P(a,b)) < 87 + 25 

as desired. 

Appendix B: Proof of Lemma [3] 



Suppose |M| = fc, i.e. M = {Mi, M2, . . . , Mfe}. We 
first prove that 



^/t(M,)/T(M,) </. 



(Bl) 



(When we write A < B for two matrices _B, we mean 
that i? — A is positive semidefinite.) Note that by the 
definition of fx {Mi), we can write it as 



MM,) ^ friM,) (g> I 



(B2) 



where friMi) is some operator on T, and / is the identity 
operator on T'^. Then Eq. (|Bip is equivalent to 



(B3) 



^/t(M,)/T(M,)</. 



i=l 



Let IV') be an arbitrary pure state on T, and p be the 
uniformly mixed state on T'^. Then by plugging 



M, ^ friM,) ^ I + g{M,) 



into 



1 = tr(|V')(V'|'^p) 

= tr[(EM/MO(|^)(VI®p)], 



(B4) 



(B5) 



and noting that tr(gT{Mi)) = for any z, we obtain 

1 = {MiE fkM^)fT{M,))m 

+tr[(X:4(M,)gT(M,))(|^)(V^| ®p)] (B6) 

> (V'I(E/;(m.)/t(m,))|v). 

1=1 

Since lip) is arbitrary, we get 



J2fkM^)fTiM,)<I 



(B7) 



i=l 



(B8) 



as desired. 

Now let TV = {TVi, 7V2, • ■ • , Nk+i}, where 

TV, -/t(M,), Vz = l,2,...,fc; 
Nk+i = EfliMi) friMi). stem (B9) 

(For a positive semidefinite matrix A with spectral 
decomposition A = Ex^^lxjixl, we define :— 
J2x Then TV is a valid measurement. More- 

over, since each friMi) acts non-trivially only on sub- 
system T, so does each TV^. Hence TV is a |T|-local mea- 
surement. 

Now, the distance between Mi and TVi satisfies 



A2(TW„TV,) < J-\\M,-N,\\l 



1 

^ll5T(M,)||| 

^imrF-\\fTiM,)rF), 



(BIO) 



for i = 1, 2, . . . , fc, and 



1 



AiMk+i,Nk+i) = ^\\Nk+,\\l 

2 - 2^ E ll/T(Af,)lll^. 
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So the distance between M and N satisfies 
A'(Af,W) = '|;'A=(Af,,iV,) 

Here we have used the fact that 

k / ^ \ 

i=l \i=l / 



Appendix C: Proof of Lemma [4] 

The proof of this lemma is similar to that of lemma |31 
Suppose \M\ = k, i.e. M = {Mi, Ma, . . . , Mfe}. We 
first show that 



or equivalently, 

fc 



(CI) 



(C2) 



Fix any A e ^d.n- Let be an arbitrary pure state in 
Wx.d, and p be the uniformly mixed state in V^. Then 
by plugging 



where 



M, = (M, + A'U + M,)geh 



= (Af,)A®/, 
M, = E |A)(A'|® (M,)a,a' 



(C3) 

(C4) 
(C5) 
(C6) 



into 



1 = tr[(|A)(A|®|V;)(V'|®p) 



SchJ 



= tr 



E m/m, (|A)(A|®|^)(V'|®p) 



Sch 



and noting that tr{g\j) = for j 7^ 0, we obtain 

1 = mtmlrnxm 



+tr[(X:M,^M,)(|A)(A|§5|V)(^|®p)] 



(C7) 



(C8) 



Since \ip) is arbitrary, we get 



E(M.)1(M.)a < /, 



(C9) 



as desired. 

Now let N = {Ni,N2, Nk+i} where 



N, = (M,)sch, Vz = l,2,...,fc; (CIO) 

(Cll) 



/-E(^^) 



Sch- 



i=l 



Then iV is a valid measurement. In addition, since each 
Mi is permutation-invariant, so is each Ni. Thus, is a 
permutation-invariant measurement . 

Now the distance between M, and iV,- satisfies 



A2(M„iV,) < -i-||M,-A^,||| 



1 

¥ ~ - 

— ||(M,-|-M,)sch||| 

— (||M,|||, + ||M,)|||,) 

— (||M,|||-||M,|||), 



(C12) 



for i = 1, 2, . . . , fc, and 



1 



A2(Mfe+i,7V,+i) < —\\Nk+i\\l 



So the distance between M and N satisfies 

k+i 



(C13) 



A2(M,iV) = EA2(M„7V,) 
i=i 



(C14) 



Here we have used the fact that 

k / fc 



Appendix D: Proof of Lemma [5] 

For convenience, let I — {i : Li > 0.16^L/k} and 
= \{v{M,)\i){Ni))\ for any i. Then by A(M,iV) > S we 
get 

52 < A2(M,7V) 

= 1 - E v^(M^M]Vi)e,. 
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It follows that 



1-6^ > E ^P{M^)p{N,)^.^ 

> (1-0.1,52) E 



(D2) 



sm.ce p{Mi),p{Ni) > (l-0.l52)Lj/L for any i e I. Hence 



(1-n^ 

0.1(52 



< (l-0.9,52)i. 



(D3) 



Meanwhile, for any i ^ I, we have Li < 0.16^ L/k, and 
there are at most k such i's, so 



E > (1 - 0.1^2)^^ 



(D4) 



Furthermore, by the fact that the arithmetic mean of 
a set of non-negative real numbers is no less than their 
geometric mean, we get 



< 



Now, by Eqs.(|D3 l) - ((D5|) and < < 1, we obtain 



(D5) 



i(xL(Af)ixL(A^))i = n^"^- 



iei 



< 



1 - 0.9<52 



2x (1-0.1<5")L 



(D6) 



1-0.1,52 

< (1 - 0.8,52)"-9^ 

< (1-0.6,52)^. 



Appendix E: Proof of Lemma [6] 

Without loss of generality we assume that 
\4'2)t---, \4'm) are linearly independent. Then the 



subspace spanned by them is m-dimensional. Let 
|(/?2), ■ • ■ , Iv'm)} be an orthonormal basis for this 
subspace, where 



i=i 



i=ij=i 

1 m 7n 

^ i=ij=i 



Hence 



■IHI^ 



m m 

= EIEA*,(0,|^)|2 
i=i 3=1 

m m 

mm /I 

< E E |A.,.l|A.,,n- 

i=lj.j'=l \^^/ 

1 ™ ™ |A,,,.|2 + |A,_^.,|2 



< 



E E 



1 771 m 

^ 25;^^^"^-"' 

< 0.1. 



(El) 



for some Xij's. Then 



i=l 

m m 

= E E A*,A.^,,(0,|,/.,-) 

= E EIA,,P + E E K3hf{<f>M 
1=1 j=i 1=1 j^j' 

mm m |\ . .12 , |\ . 12 -i 

> EEIA..P-E E V 



(E2) 



(E3) 
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